Pauli spin blockade and the ultrasmall magnetic field effect 
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Based on the spin-blockade model for organic magnetoresistance we present an analytic expression 
for the polaron-bipolaron transition rate, taking into account the effective nuclear fields on the sites. 
We reveal the physics producing qualitatively different magnetoconductance line shapes as well as 
the ultrasmall magnetic field effect, and we study the role of the ratio between the intersite hopping 
rate and the typical magnitude of the nuclear fields. Our findings are in agreement with recent 
experiments and numerical simulations. 



The discovery some ten years ago of spin injection in 
organic semiconductors [T] and a giant magnetoresistance 
in organic spin valves |H [3] triggered the birth of the 
thriving field of organic spintronics 4 , which offers in- 
teresting new physics and the potential of industrial ap- 
plications. An exciting phenomenon in this field is a large 
(up to 20%) magnetoresistance observed in different or- 
ganic materials [5HZ]) usually at small magnetic fields (1- 
10 mT) but sometimes at larger fields (10-100 mT), and 
persisting up to room temperature. Since its discovery in 
2004, different explanations for this organic magnetore- 
sistance (OMAR) have been proposed: For bipolar de- 
vices it was suggested that spin-dependent electron-hole 
recombination and dissociation rates could be responsible 
[3] , whereas a model based on nuclear-field-mediated 
bipolaron formation could explain OMAR in both bipo- 
lar and unipolar devices [Hflfi"2] . 

More recently an organic magnetoresistivc effect on an 
even smaller field scale (0.1-1 mT) has been observed in 
unipolar as well as bipolar devices [13]. This ultrasmall 
magnetic field effect (USMFE) is manifested by a sign 
reversal of the magnetoconductance (MC) at very small 
fields, creating two small peaks(dips) around zero field 
for devices with a negative(positive) MC. Experimental 
results seem to indicate that the typical field magnitude 
on which the USMFE is observed scales with the width of 
the MC curve when different materials are investigated 
[13] . An explanation for the effect was suggested in terms 
of enhanced singlet-triplet mixing close to the crossings 
of the hyperfine sublevels of pairs of charge carriers (po- 
larons) coupled to single nuclear spins [14] • This expla- 
nation is still under debate, mainly because it is expected 
that a single polaron in reality couples to many nuclear 
spins [15H17) . Numerical simulations based on a semi- 
classical model (where the coupling to an ensemble of 
nuclear spins is treated as an effective magnetic field) 
also reproduce the USMFE [IS] and thus invite to seek 
for an explanation along semiclassical lines. 

Here, we study the OMAR line shape as it naturally 
emerges from the spin blockade model of Ref. [10] . We 
present an analytic expression for the charge current 



through a polaron-bipolaron link for a given realization 
of the nuclear fields. Our results reproduce the USMFE 
and the different line widths as observed in experiment 
and in numerical calculations based on the same semiclas- 
sical approach [18] . and from our analytic insight we can 
identify the underlying physical mechanisms. We note 
that many interesting aspects of spin-blockade physics 
have already been investigated in the seemingly foreign 
field of spin qubits hosted in semiconductor quantum 
dots [HI [20] i where spin blockade is commonly used as a 
tool for single-qubit readout [21] [22] . Indeed, the physics 
of the polaron spin blockade model for OMAR is very 
similar to that governing the electron transport through 
a double quantum dot in the spin-blockade regime |20) . 
Our investigation thus builds on the theoretical frame- 
work of Ref. [20], and our explanation of the USMFE 
relies on a subtlety which was not addressed in Ref. [20] . 

Let us first briefly review the bipolaron model for 
OMAR presented in Ref. [10]. Electric current flows 
through the organic material as polarons hop between 
different localized molecular sites. Typically, the sites 
participating in transport do not form a regular lattice 
and all have a random energy offset with a distribution 
width it of 0.1-0.2 eV [ID]. Sites with a relatively large 
negative energy offset are likely to trap a polaron for a 
long time, but since the on-site polaron-polaron repulsion 
is typically of the same order of magnitude as cr, such 
occupied sites can often still take part in transport by 
temporarily hosting a pair of polarons, i.e. a bipolaron. 

Due to a relatively large orbital level spacing, most en- 
ergetically accessible bipolaron states are spin singlets. 
This makes the polaron-bipolaron transition spin selec- 
tive, ultimately leading to OMAR. The mechanism can 
be understood from Fig. [l] where we focus on a single 
polaron-bipolaron transition. We assume that the spins 
of the two encountering polarons are random and for sim- 
plicity we describe the problem in the basis of spin eigen- 
states quantized along the direction of the local magnetic 
fields B^ij. Two possible initial spin states are depicted: 
(i) the left spin antiparallel and the right spin parallel to 
the local field, and (ii) both spins parallel. In the absence 



2 




FIG. 1. (color online) When a site has a relatively large neg- 
ative energy offset and already contains a polaron (the right 
site in the pictures), charge transport through this site relies 
on the formation of a bipolaron (blue arrows). Large on-site 
exchange effects dictate this bipolaron to be a spin-singlet, 
which leads to spin-blockade physics. 



of an external field, the magnetic fields at the two sites 
are the local random effective nuclear fields [Fig. [ija)]. 
Generally all initial states can then transition to a spin- 
singlet bipolaron and current runs through the system, 
ff, on the other hand, a magnetic field much larger than 
the typical nuclear fields is applied, then and B# are 
(almost) parallel [Fig. [ljb)]. fn this case, situation (ii) is 
a spin triplet, out of which a bipolaron cannot be formed: 
the current is blocked. 

We thus see that a simple two-site picture is able to 
explain the essentials of OMAR. Of course, in experiment 
there are many possible paths for charge carriers through 
the material and not all of them contain bipolaron sites. 
The visibility of all effects of spin blockade will thus be 
reduced, but the characteristic features survive [TU] . 

In this work, we will focus on the physics of a single 
polaron-bipolaron transition and its MC line shape. To 
describe the transition, we use five states: the four pos- 
sible initial spin states of the polaron pair (both sites 
hosting one polaron), one spin-singlet state \S) and three 
spin-triplet states |T ) and \T±), and the spin-singlet 
bipolaron state |Sb). The Hamiltonian we use to describe 
the coherent dynamics of these states reads [5D] 
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written in the basis {\T + ) , \T ) , |T_) , \S) , \S h )}. This 
Hamiltonian includes a coupling energy t between the 
two singlets (which enables polaron hopping) and the rel- 
ative energy offset (detuning) A of the bipolaron state, 
typically A ~ a. The effect of the local magnetic fields 
is expressed in terms of the sum and difference fields 
B s = |(Bx, + B R ) and B Q = |(Bi - B R ), and we use 
the notation B , > = , ± iB y , ,). Note that we 

s{a) v/2 v s{a) s(a)' 

have set gps = 1 for convenience. 



As pointed out in Ref . [5D] , we can deduce already from 
Eq. ([!]) that there exist in the space of (B^, B^) so-called 
"stopping points" where the current is blocked. To see 
this, we take the spin quantization axis to point along B s , 
which amounts to setting Bf — > in Q. Then we find 
that current vanishes when B a || B s or B a _L B s , since 
at these points one or more of the triplet states are not 
coupled to \S). The sum and difference fields B s a both 
contain a contribution from the effective nuclear fields 
Ki /j on the two sites, whereas the external field B cxt 
only adds to the sum field: B s = K s + B cyi tz and B Q = 
K a . For a given random realization of I£l,r one can thus 
always find a field B oxt for which B a _L B s , and a sweep 
of -Bext for a fixed K-l,r will always exhibit a stopping 
point where the current vanishes. The position of this 
stopping point is determined by the relative orientation 
of K s and K a and is thus random. In an experiment one 
usually sweeps B cxt so slowly that at each measurement 
many configurations of the fields Ki^ are probed. As a 
result the stopping points are averaged out and one finds 
a smooth MC curve [20] . 

However, this is not the full story. A subtlety, not 
discussed in Ref. [20], is that there exists one more stop- 
ping point [23 : When B s = the triplet subspace in 
the Hamiltonian is degenerate and Eq. ([I]) can be equiv- 
alently written in terms of one coupled triplet state 

._ . _ -B- \T+) + Bl \T ) + B+ |T_) 



and two orthogonal triplet states \T\) and IT2) which 
have (Ti^\H\S) = and are thus blocked. Why would 
we bother? We argued above that stopping points oc- 
cur at random positions, leaving no trace after averaging 
over Ki_j;. This new stopping point however, is fun- 
damentally different from the ones discussed above: It 
suppresses current close to where B oxt = — K s , which is 
always in the vicinity of -B C xt = 0. It is therefore possible 
that after averaging over K^ jj this new stopping point 
leaves a trace in the MC curve: A small dip around zero 
field, like the USMFE. 

Let us now explicitly calculate the current as governed 
by this polaron-bipolaron transition. To describe charge 
transport, we write a time-evolution equation for the 5x5 
density matrix of the system. To the coherent evolution 
dictated by H we add incoherent rates describing disso- 
ciation of the bipolaron to the environment and hopping 
of a new polaron onto the empty left site. This yields 
(where we have set h to 1) 

^ - -i[H, p] - ^{P h ,p} + Jp b ,b(l - (2) 

where T is the rate of bipolaron dissociation to the envi- 
ronment, P b = \Sb) (Sb| is the projection operator onto 
the bipolaron state, and 1 is the identity matrix. In 
writing so we assumed for simplicity that refilling of the 
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left site takes place immediately after dissociation of the 
bipolaron. If this is not the case, the prefactor for the cur- 
rent changes but the MC characteristics stay the same. 

We add the normalization condition Tr[p] = 1 to the 
set of equations and then solve dtp^ c °^ = to find the 
stationary density matrix. The charge current is then 
given by I = eTp/j^ and can be found explicitly. We as- 
sume for convenience that r 3> t, B s , B a , A is the largest 
energy scale in the problem [T2], and then find 
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(3) 



in terms of x = B s /B a . Here Y s = t 2 /Y is the singlet- 
singlet hopping rate from the left to the right site and 4> 
is the angle between B s and B a [24]. We also used 
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2 cos 2(f>. 



(4) 



We see that all stopping points predicted above are in- 
deed reflected in At B s = we have x = which 
yields I — 0, and B Q |j B s or B a _L B s corresponds to 
<j> = 0, 7T or <f> = 7r/2 respectively, both also giving 7 = 0. 

Equation |3]) is the most important analytic result of 
our work. It gives the current for one single realization 
of Ki, Kr, and i? cx t- The MC measured in experiment 
is found by averaging ^ over the random nuclear fields. 
In contrast to the analytic results presented in [2U], our 
result is valid for arbitrary Y s and not only for limiting 
cases. One word of caution is required here concerning 
the interpretation of ([3]): If one wants to plot I(B cxt ) for 
a single realization of r, one should not only use B s = 
| K s + .Bext^l m ^ but also implement the dependence 
of <p on B cxt implied by cos cf) = (B s ■ B a )/B s B a . 

Let us now investigate Eq. ^ and see what we can 
infer about the line shape of the predicted MC curve. 
We always have a > — 2, which ensures that I > ev- 
erywhere. The current vanishes for x — or x — > oo, 
and in the range x S [0, oo] we have a single maximum 
at x = 1 where the current is / max = 4er s sin 2 (f>/(a + 2). 
In Fig. [2]ja) we plot the expression given in Eq. ^ for 
different a. The FWHMs w_ of the dip around x = and 
w + of the overall peak structure (as indicated in the plot 
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for a = 4) are found to be w± = 2+ |a± y 3 + 2a 

We see from Eq. Q that an important parameter is 
Y s /K, the ratio of the intersite hopping rate and the 
typical magnitude of the nuclear fields K, typically ~ 
0.1 /xeV [T5] • We will thus now investigate the cases of 
small and large Y s /K. 

In the limit of Y s /K <C 1 we can write 
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= Y s (n L xn R ) 2 , (5) 



where we used the unit vectors n^.R = Blm/Bl^. As it 
should, this result coincides with that of Ref. 201 in the 
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FIG. 2. (color online) (a) The current given by (J3j) as a 
function of x = B s /B a for fixed B a and 0, evaluated for 
different parameters a. (b,c) The averaged MC. (b) Blue 
trace: Y s /K <^ 1. The field B xt is plotted in units of 
K = {Kl, R ) 1/2 . The peak of this curve is flat [20]. Green 
trace: Y s /K = f . (c) T s /K = 50. Now B cxt is plotted in 
units of T s . (inset) The range where _B cxt ~ K. 



same limit: There are no intersite exchange effects and 
the situation is exactly like the picture of Fig. [T] where 
the current only depends on the relative orientation of 
Bl and B^. As was shown in [20 , Eq. ^ can be av- 
eraged analytically over random K.l.r taken from a nor- 
mal distribution, yielding a MC curve with a flat peak at 
_B cxt = 0, a maximum of (I) max ~ eY s , and a line width 
of ~ K. Indeed, for all a € [—2, 2] we find that w+ ~ 1, 
so for any <fi the current is suppressed when x ^ 1. In 
Fig.[2]jb) (blue trace) we plot the resulting MC line shape, 
where we defined MC(B cxt ) = [I(B cxt ) - 1(0)}/ 1(0). 

In the opposite limit of Y s /K ^> 1 we have a ps 
(r s /B a ) 2 (3 + cos -2 <f>) ^> 1. We can already see from the 
properties of Eq. (I3J) that in this case (I) max ~ eK 2 /Y s , 
and that w+ ~ aJi 2 ~ r s /i4T implies a MC line width 
of ~ r s . Indeed, T s sets the level broadening of \S) and 
as long as B s < r s generally all three triplet states can 
efficiently transition to \S) with the coupling provided by 
B Q . The width of the dip around the stopping point at 
x = is w_ a -1 / 2 ~ K/Y s in terms of x, or ~ K 2 /Y s 
in terms of B s . This energy scale can also be understood: 
If B s = the decay rate of \T m ) is Y t ~ K 2 /Y s , which is 
the only energy relevant in the triplet subspace. When 
B s ^ Y t the decay of the other two triplet states becomes 
comparable to Y t and the blockade is lifted. 

In this limit of large Y s /K the current cannot be av- 
eraged analytically over the nuclear fields, and one has 
to evaluate the integrals over the distribution of Hl.r 
numerically. Fig. [2jc) shows a plot of the MC integrated 
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| sweep of Z? cxt 




FIG. 3. When _B ex t is swept for a given realization of K_l,_r, 
the field B s (blue arrow) follows a trace like the blue dashed 
line. We indicated with green arrows K s as well as , which 
equals the minimum value of B a . 



over normal distributions for all six components of 
and Kr. For all components we used a standard devia- 
tion of K/y3 and we have set T s /K = 50. The resulting 
line shape is Lorentzian since it is determined by the level 
broadening of \S). Close to zero field, where B ext ~ K, 
we find a very faint USMFE, as shown in the inset. When 
we set T s /K even larger we find that the visibility of this 
USMFE is suppressed further, ultimately reaching zero. 

We can understand this USMFE from the expression 
for the current given in Eq. ^, For a given realization 
of ~Kl,r, the current trace I(B cxt ) "misses" the zero-field 
stopping point by K^- = \J (Ks) 2 + (Ks) 2 , as illustrated 
in Fig. m Some realizations have > B a so that the 
current trace has a single maximum [see Fig.^a)]. Other 
realizations have < B a and the current exhibits a dip 
at small fields, the position of the dip at B cxt = —K*. For 
large T 8 /K the dip around the zero-field stopping point 
becomes narrow, of the order ~ K 2 /T s <C K , and only 
the very few curves of I(B cxt ) with K~- < K 2 /T s have an 
appreciable dip. This still can produce a faint dip in the 
averaged current. However, the position of each single- 
realization narrow dip is at B ext = so averaging 
over K z s makes the averaged dip even less pronounced 
and results in a dip width of ~ K. 

The regime to look for a pronounced USMFE is thus 
at intermediate T s /K ~ 1. In Fig. [2jb) (green trace) we 
plot the averaged MC for T s /K = 3/2 and we see indeed 
a strong USMFE, its visibility being ~ 5%. This regime 
is optimal for the USMFE since here the width of the 
zero-field dip is still ~ K but the symmetric situation 
where the current only depends on the angle between 
wl and riR is significantly perturbed. In other words, 
at T s /K — > the overall MC line width is minimal and 
~ K. The two USMFE "bumps" are still there but are 
split by the same energy scale ~ K and thus appear just 
left and right of the top of the MC curve. In the limit 
oiT s /K = the bumps and the underlying MC curve 
have exactly compatible shape line shapes, together re- 
sulting in the characteristic flat peak. If one moves away 
from T s /K — the underlying MC line shape becomes 



broader, which makes the USMFE bumps more visible. 
However, as soon as T s /K becomes too large, one enters 
the regime discussed above, where the USMFE disap- 
pears again. The optimal regime is thus at T s /K <~ 1, 
in agreement with the results presented in Fig. [2]jb,c) as 
well as with previously obtained numerical results |12j . 

To summarize, we studied the two-site spin-blockade 
model for OMAR and derived an analytic expression for 
the polaron-bipolaron transition rate, taking into account 
the local nuclear fields on the two sites. We showed how 
our result reproduces different MC line widths: ~ K for 
slow and ~ V s for fast intersite hopping. We also pro- 
vided an explanation of the USMFE in terms of a persis- 
tent spin blockade at the special point where the average 
effective field vanishes B s — 0. The USMFE as predicted 
here always takes place on the field scale ~ K, and we 
explained why it is expected to be most pronounced in 
the regime where V s ~ K . In this regime thus both the 
scale of the USMFE and the MC line width are set by 
K, the latter however being slightly larger. This relation 
between the two scales is consistent with experimental 
observations |13j and numerical simulations |12i . 

As a side remark we note here that a close inspection 
of the experimental data presented in Ref. [20] (the cur- 
rent through a double quantum dot) also seems to reveal 
a faint USMFE. The data were fitted to the flat-peak 
curve since the system was assumed to be in the inelas- 
tic tunneling regime. In reality the coupling is however 
never perfectly inelastic, and a faint trace of the USMFE 
could be left. Due to its tunability, a double quantum 
dot might in fact be the best system to experimentally 
explore USMFE in more detail. 

We acknowledge helpful feedback from M. S. Rudner 
and P. A. Bobbert. 
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